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In order to clarify the physics of the gating of solids by ionic liquids (ILs) we have gated lightly 
doped p-Si, which is so well studied that it can be called the ’’hydrogen atom of solid state physics” 
and can be used as a test bed for ionic liquids. We explore the case where the concentration of 
induced holes at the Si surface is below 10 12 cm~ 2 , hundreds of times smaller than record values. 
We find that in this case an excess negative ion binds a hole on the interface between the IL and 
Si becoming a surface acceptor. We study the surface conductance of holes hopping between such 
nearest neighbor acceptors. Analyzing the acceptor concentration dependence of this conductivity, 
we find that the localization length of a hole is in reasonable agreement with our direct variational 
calculation of its binding energy. The observed hopping conductivity resembles that of well studied 
Na + implanted Si MOSFETs. 


I. INTRODUCTION 


In the last decade room temperature ionic liquids rev¬ 
olutionized the way in which weakly conducting or insu¬ 
lating materials are gated in order to make them metallic 
or superconducting. By using an ionic liquid (IL) one can 
achieve surface carrier densities of n ~ 10 15 errD 2 within 
the electrochemical window of the electrolyte, which is of 
the order of a few volts P 

Ionic liquids are molten salts made of negative and 
positive ions with relatively large radii. As a result, the 
Coulomb attraction energy of two adjacent opposite sign 
ions, Ec = e 2 / kD, is small enough so that at room tem¬ 
perature the salt stays molten. Here D ~ 1 nm is the sum 
of radii of positive and negative ions and k ~ 3 is the di¬ 
electric constant of the IL. A negative voltage, applied 
between a metal coil submerged in the IL and the sam¬ 
ple, causes negatively charged ions (anions) to migrate 
to the sample surface. At the sample surface an electric 
double layer is formed, which can be viewed as a capac¬ 
itor with nm scale separation between electrodes. The 
crucial advantage of IL gating is that the IL-Si junction 
is insulating. 

So far most of the research has dealt with gat¬ 
ing of novel materials, such as rubrene 3 , underdoped 
YBCO S, polymersP, nanocrystal arrayiP, nanotubes, 
and graphene^, and has been aimed at the record induced 
carrier densities. 

In this paper we focus on the basics of gating with ILs. 
For this purpose we choose to gate a lightly doped single 
crystal p-type Si wafer. Silicon is so well studied that it 
can be called “the hydrogen atom of solid state physics” 
and, therefore, used as a “test bed” for ILs. At large sur¬ 
face densities the IL charge layer is usually regarded as an 
uniform sheet. In this paper we consider relatively small 
negative gate voltages, when the concentration of surface 
holes is in the range of 10 11 cm -2 < n < 10 12 cm -2 . In 
this case each excess ion driven to the surface should be 
treated as a discrete charge individually responsible for 
creation of a carrier in a sample. 
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FIG. 1. Surface bound states of carriers in Si (a) gated with 
IL and (b) gated by Na + implantation in the oxide. In a) 
the red region represents the wave function of the hole. The 
negative charge acceptor cluster of 3 ions is shown by blue in 
a) at the distance d from the Si surface. In b) the electron 
cloud is shown by blue, a is the decay length of the surface 
bound state wave function in Si. The spherical coordinates 
r, 0 are used below for the variational calculation of the hole 
binding energy. 


In the absence of an applied gate voltage, the IL is a 
neutral and strongly correlated liquid. When Ec fceT 
the arrangement of positive and negative ions resembles 
that of NaCl. When the gate is negatively biased, ad¬ 
ditional anions are driven to the sample surface. This 
results in discrete excess negative ions imposed on the 
neutral background with NaCl-like order. An example 
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of such a net negative group of ions at the surface of Si 
is shown on Fig. Below we present evidence that 

suggests the excess anion binds a hole residing in Si (see 
the hole wave function cloud in Fig. EM* One can no¬ 
tice that such a IL surface acceptor is similar to the Na + 
surface donor in Si MOSFETs as shown in Fig. 

In the Na + experiments, the Na + ions were implanted 
on the top of layer of Si 02 and driven by a positive gate 
voltage to the Si surface. Each Na + ion binds an electron 
near the Si surface, becoming a surface donor. The sur¬ 
face conductivity as a function of temperature was shown 
to follow 
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where £j is the activation energy and cq the conductivity 
prefactorP The first term, which represents the activa¬ 
tion of electrons to the conduction band with the activa¬ 
tion energy e\ equal to the ionization energy of the donor, 
dominates at high temperatures. The second term, ob¬ 
served at low temperatures, describes the nearest neigh¬ 
bor hopping between the Na + donor states at the Si sur¬ 
face. (Some compensation of Na + donors by acceptors 
creates empty donors and allows nearest neighbor hop¬ 
ping to proceed.) 

In this work we show that IL gating of p-type Si leads to 
a similar transport phenomenon. We observe the nearest 
neighbor hopping with the activation energy £3 changing 
with the gate voltage, while the first term of Eq. ([T]) 
is shorted by the bulk conductance. This confirms the 
existence of discrete and sparse IL acceptors at the Si 
surface shown in Fig. HI 1 - 

The remainder of this paper is organized as follows. In 
the following Section [II] we discuss our experimental set 
up. In Section m we present our experimental results 
and their analysis. In Section IV we discuss our model 
of IL gating of p-Si with a Pt electrode and explain the 
origin of the negative threshold voltage. In Section [V] we 
theoretically evaluate parameters of the surface acceptor 
and compare them with our experimental results. We 
conclude in Section IVTl 


II. EXPERIMENTAL DETAILS 

In our experiment we made ohmic contacts to a boron 
doped Si wafer after annealing it in a mixture of N 2 and 
5% H 2 in order to reduce the number of dangling bonds 
at the surface. A glass cylinder attached to the Si surface 
was used to confine the ionic liquid and control the gated 
area between contacts. The gate electrode consisted of a 
Pt metal coil suspended in the ionic liquid DEME-TFSI. 
The resistivity of the wafer is 1 — 5 flcm which with the 
typical hole mobility in Si of 400 cm 2 /Vs corresponds 
to a boron concentration of 10 16 cm . This concen¬ 
tration is smaller than 3D concentration of 10 1 ' cm -3 
( 10 11 cm - 2 / 10 -6 cm) induced by the gate in the surface 


layer which has a width of 10 -6 cm (see below). Addi¬ 
tional experimental details can be found in the Supple¬ 
mental Material of Ref. m 

Other experiments have reported the possibl e rol e of 
electrochemical reactions when gating with an IL 111 . To 
avoid these reactions we modulated the carrier density 
below room temperature at T = 230K, where the IL ions 
are still mobile. A thin oxide was used to protect Si from 
further oxidation and in ad dition , to passivate the trap¬ 
ping sites on the Si surfac^mi_ This native oxide has 
a thickness of 10 — 15 A as measured by ellipsometry (a 
thicker oxide would reduce the capacitance of the tran¬ 
sistor) . 

The 4-terminal sheet resistance Rs of the silicon chan¬ 
nel was measured by the van der Pauw method. In our 
previous stud ;/ 10 we showed that the temperature depen¬ 
dence of the resistance at each gate voltage was repro¬ 
ducible during the cooling to 2 K and warming up to 230 
K. 

To calculate the surface concentration one can, in prin¬ 
cipal, use the capacitance of the device. However, in such 
a structure the total capacitance depends on both the 
geometrical and quantum capacitancdA 4 . The last one 
can be calculated for a degenerate gas with high hole 
concentratiorP^l or for small concentration^ 15 -! But our 
situation is one of intermediate concentration when the 
distance between holes is comparable with the Bohr ra¬ 
dius of a hole. In this case one cannot use these formulas. 
Doing so would lead to uncertainties in the capacitance 
value. Therefore, to calculate the surface hole concentra¬ 
tion we rely on integration of the charging current. 

Just like in the Na + experiments we do not observe 
the surface conduction until a certain sample dependent 
threshold gate voltage Vt = —0.7 V is reached^. At 
a gate voltages Vg lower than Vr an additional anion 
brings and binds a free hole from the bulk of the Si sam¬ 
ple, becoming a shallow acceptor. Therefore, the two- 
dimensional hole density, n, at Vg < Vt was calculated 
by integrating the current during the charging process 
starting from Vg = Vt- The dependence of the inte¬ 
grated carrier de nsity n on the applied gate voltage vjsn 
is shown on Fig. [ 2 p£l 

From the slope of the dependence of n on Vq one can 
find the capacitance of the device 2 ^F/cm 2 . Assuming 
that dielectric constant of the IL and the SiC >2 is the same 
kso = 4 the related distance between negative ions in the 
IL and the holes in Si is 17 A . This agrees reasonably 
with the sum of the size of the ion and the width of the 
SiC >2 layer. 

Since our device has a large area ~ 6 mm 2 and we drive 
relatively large current ~ 50 nA the noise contribution 
to the accumulated charge is negligible. 


III. RESULTS 

Plotted as a function of inverse temperature, the sheet 
resistance Rs is shown in a logarithmic scale in Fig. [3] 
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FIG. 2. Integrated carrier density n as a function of applied 
gate voltage Vg- 

at different Vg corresponding to four different concen¬ 
trations of holes. Above 20 K, Rg plotted in this way 
follows a straight line independent of Vg 

Rs = R » exp (if?)’ (2) 

where e& = 43 meV and Rq is a prefactor. The ionization 
energy of a Boron acceptor in Si is known to be 45 meV, 
so that it is natural to assume that this line corresponds 
to the bulk conductance. Below 20 K, Rg continues to 
follow activated temperature behavior, but with a Vg 
dependent activation energy. 

It is natural to interpret this low temperature part of 
Fig. [3] as the nearest neighbor hopping conductivity of 
holes between IL acceptors, which can be described by 
the second term of Eq. 0. 

In our p-type Si samples 1/Rq a i; therefore, the 
(7 1 part of the surface conductivity of Eq. ([!]) is shorted 
by the bulk conductance. This shorting does not occurs 
in Na + doped samples because they are in the inversion 
layer regime, where the conducting path between source 
and drain electrode is insulated from the bulk of the sam¬ 
ple by the hole depletion layer, while in our case we deal 
with an accumulation layer. 

From each plot of Fig. [3] we extracted the prefactor 
<73 and the activation energy 63 which are shown in Fig. 
[4] and Fig. [5] respectively, as functions of the carrier 
density n. For comparison we added the corresponding 
data for the nearest neighbor hopping conduction of Na + 
implanted Si MOSFETs as summarized in Table III of 
Ref. [9} Remarkably, the two experiments cover the same 
range of carrier concentrations and show similar hopping 
parameters. 

The dependence <73 (n) can be compared with the theo¬ 
retical prediction based on the percolation approadPUfUl 
as explained in Chapters 5 and 6 of Ref. [2jJ, The two- 


FIG. 3. Natural logarithm of the sheet resistance R s (in 
Ohms) plotted as a function of inverse temperature. Bulk 
conduction is independent of the gate voltage. When the ions 
create a surface channel, nearest neighbor hopping is observed 
at low temperatures. The four curves correspond to different 
surface concentrations n, in units 10 11 cm -2 . Two carrier 
densities are omitted for clairity. 



FIG. 4. Prefactor of nearest neighbor hopping conductivity 
vs. the inverse square root of the carrier density of surface 
carriers induced by IL gating (this work) and by Na + implan¬ 
tation (Ref. [9j). Straight line fits are drawn to evaluate the 
decay length a of the localized carrier wave function on IL 
acceptors and Na + donors as explained in the text. The size 
of the symbols reflects experimental uncertainty. 

dimensional analog of Eq. (6.1.9) with the help of Table 
5.3 gives 

<73 = 030 exp[—2.4/Vna 2 ], (3) 

where (730 is a conductivity prefactor. In Fig. [4] we 
compare our data with Eq. [3] A reasonably good fit 
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of the data is provided by the two straight lines corre¬ 
sponding to this equation. For the decay lengths we get 
a = 59 ± 10A for Na and a = 83 ± 20A for IL. The rela¬ 
tively large error in a for the IL case comes from the ac¬ 
curacy of our carrier density which was determined from 
the uncertainty in Vp to be An = 6x 10 10 cm -2 . 


IV. MODEL OF IONIC LIQUID GATING 

Let us present our model of IL gating in greater de¬ 
tail than in the introduction. Let us discuss what hap¬ 
pens when a Pt electrode gates p-Si crystal with help 
of an IL but still at Vq = 0. The work function of 
Pt, ippt — 5.7 eV, while the work function of p-Si, 
<PSi = X + E g = 5.1 eV, where x = 4 eV is the elec¬ 
tron affinity and E g = 1.1 eV is the band gap of Si. In 
equilibrium the electrochemical potentials of Pt and p-Si 
have to be equal. This means that electrons are trans- 
fered from p-Si to Pt. This results in a positively charged 
p-Si and a negatively charged Pt. The resulting electric 
field between them is compensated by the IL polariza¬ 
tion, so that anions move to p-Si and cations move to the 
Pt surface. The voltage difference ppt — <psi — 0.6 eV 
drops at the capacitance of the IL/Si surface, because 
the area of the Pt electrode is much larger than that of 
the Si one. Thus, there is an excess number of anions on 
the Si surface even in an equilibrium. This means that 
ideal contact Pt/IL/Si accumulates a layer of holes in Si. 
In such a situation, the Si wafer surface would be a good 
p-type conductor. 

In fact the situation is different. An ungated p-type Si 
wafer is covered by a thin oxide layer. It is known that 
the oxide-silicon interface contains a substantial concen¬ 
tration of oxygen vacancies, which play the role of sur¬ 
face donors in p-Si 24 . This results in an inversion layer of 
electrons on the ungated surface of p-Si. In other words, 
the Si bands are bent downwards by the value of the Si 
gap Eg creating a pocket for electrons while holes are 
depleted in a wide layer near the surface. 

Let us see what happens when such a real surface is 
gated by an IL with a Pt electrode. In this case, p-Si 
donates some of its electrons to Pt and becomes posi¬ 
tive. As result of the IL polarization an equal number of 
anions move close to the Si surface forming a capacitor 
with positive donors. Actually at the temperature 230 K 
where gating is done the Coulomb attraction energy of 
an anion to a positive donor is substantially larger than 
kpT so that most of the anions form compact pairs with 
empty positively charged donors. The remaining donors 
still create a downward bending of the conduction and 
valence bands and the depletion layer of holes, but this 
bending is smaller than E g and does not create an in¬ 
version layer near the surface. Thus, at V/; = 0 the Si 
surface is an insulator. 

To make it conducting we apply the negative Vq = 
Vp — —0.7 eV. This threshold voltage brings enough 
anions to the Si surface to neutralize all the remaining 


donors. When all the donors are paired with anions, the 
electric field applied to Si and the downward bending of 
the Si the bands near the Si surface vanish, making the 
bands flat. At Vq < Vp excessive anions bend the bands 
near the surface upwards and form a hole accumulation 
layer. Each excess anion binds a Si hole forming a surface 
acceptor discussed in the Introduction and shown in Fig. 
|T[ These acceptors overlap weakly with each other. At 
the same time their effective three-dimensional density 
easily exceeds the density of the bulk B impurities so 
that at low temperatures the nearest neighbor hopping 
of holes between surface acceptors dominates the sample 
conductance (see Fig. [3]). 

Recall that surface acceptors are situated on a back¬ 
ground of many compact pairs of donors and anions. 
These pairs are compact dipoles and create only a rela¬ 
tively weak random potential. In the following discussion 
of nearest hopping transport of holes on surface acceptors 
we ignore the potential of typical pairs. However, a small 
fraction of these pairs plays an important role in nearest 
neighbor hopping. These are pairs, which are thermally 
dissociated, i. e. where donors have their partner anion 
lost to the bulk of the IL. Such ionized donors play the 
role of a compensating impurity for the nearest neighbor 
hopping. They are responsible for the small number of 
empty surface acceptors, which allow holes to move along 
chains of surface acceptors (see Chapter 3 of Ref. 19). 


V. BINDING ENERGY AND THE HOLE 
DECAY LENGTH OF THE SURFACE 
ACCEPTOR 

In this section we theoretically calculate the binding 
energy and the decay length a for an isolated IL acceptor 
shown in Fig. d}i. Let us first review the acceptor the¬ 
ory in the bulk p-Si. It is known that the valence band 
is doubly degenerate and has two types of holes, light 
and heavy, with masses toj = 0.16 m e , to/j = 0.46 m e 
respectively, here m e is the electron mass. It is knowrPSl 
that for the above ratio mh/mi ~ 2.9 the ground state 
energy of the acceptor is E 0 = me 4 / [2K 2 Si H 2 ), where 
m = 0.7 TOft, = 0.3 m e and K$i = 12 is the dielectric 
constant of Si. If we assume that a bulk acceptor can be 
treated as a hydrogen atom with a hole mass to, then the 
wave function of the first excited p-state vanishes at the 
acceptor location and has a binding energy E 0 / 4—. 

Let us now switch to the IL acceptor. If the distance 
between the acceptor nucleus and Si surface, d , is so small 
that d <C a, one can assume, in the zero approxima¬ 
tion, that the hole is also located just near the inter¬ 
face. Then the interaction between the hole and accep¬ 
tor can be described by an effective dielectric constant 
k = (nsi + kso)/2 = 8, where kso = 4 is the dielectric 
constant of Si02- (For simplicity we make the approxi¬ 
mation that dielectric constant of the IL is that of SiC> 2 .) 
The interface between Si02 and Si is treated as an infinite 
barrier, so that the wave function of the hole vanishes at 


5 


the interface. The ground state for a hole in the half 
space is the first excited state for a hole in bulk with a 
corrected dielectric constant k. i.e. e 3 = Ry/A = 16 meV, 
where Ry = me 4 /(2n 2 h 2 ) = 63 meV. At the same time, 
the decay length a is determined by the light mass 


h 

y/2rrnei 


( 4 ) 


For ei = 16 meV we get a = 38A. 

Above we did not take into account that a hole inter¬ 
acts with its image charge and that d is not zero. To do 
so, we use a variational approach with a probe hole wave 
function 


W = 


2 


(r cos 9 — d) 
r W n (d + 2r 0 ) 


exp 



( 5 ) 


Here we use spherical coordinates shown in Fig. [l}i, r > d 
is the distance between the acceptor and the hole, 9 is 
the polar angle and the distance r$ minimizes the energy: 




'FU'FdV. 


( 6 ) 


Here the integration is over the Si half space, 9 ranges 
from 0 to 7 t/ 2 and r ranges from d/cos9 to oo. The 
Hamiltonian for the problem is 


H = - 



2e 2 

(«Si + n S o)r 

e 2 _ nsi ~ «so 

4(r cos 9 - d) Ksi(nsi + kso) ’ 


( 7 ) 


Here the first term represents the kinetic energy of a hole 
in Si. The second term is the potential energy of interac¬ 
tion of the hole with the acceptor . The third term is the 
interaction between the hole and its image in SiCU- As 
a rough estimate, we choose d as the sum of the width 
of SiC >2 layer, 10A and the radius of the negative ions in 
the IL (TFSI) 4.4A to be d ~ 15A. In the result, for 
d = 15A we get that e\ = 8 meV and a = 53A. This 
decay length for a single acceptor can be considered as 
a reasonable estimate from below to the experimental 
a = 83A. The concentrations of carriers which we work 
with are far from the case of the light doping when the 
separation between acceptors is much larger than a. The 
wave function overlap of neighboring acceptors leads to 
an additional increase of the experimental value of a. 

The fact that the overlap between acceptor states is 
substantial is also supported by Fig. [5j Here the straight 
line follows the equation: 


e 3 = 0.98—n 1/2 , (8) 

K 

which we derived for a two-dimensional lightly doped and 
weakly compensated semiconductor in a manner similar 



FIG. 5. Activation energy of hopping conductivity of surface 
carriers induced by IL gating and Na + implantation (Ref. [5} 
as a function of the square root of the carrier density. The 
sizes of symbols reflect experimental uncertainty. The solid 
line following Eq. | 8 | represents the predicted dependence of 
€3 on n at small n . 


to the derivation for three-dimensional case in Ref. 123 
(see also Chapter 3 in Ref. (23j) - Equation (|sj) takes into 
account only classical shifts of levels of holes induced by 
the Coulomb potential of rare free of IL anions surface 
donors (see end of previous section IVI and ignores the 
overlap of hole wave functions. We see in Fig. [5] that 
in contrast to Eq. [8] the experimental e 3 values decrease 
with the growing n. This happens because our concen¬ 
trations are too large . At even larger concentrations 
n > 10 12 cm -2 the activation energy e 3 vanishes and the 
system goes through the metal-insulator transition stud¬ 
ied in the Ref. m In three-dimensional semiconductors 
a larger range of concentrations was studied, which dis¬ 
played both regimes of growing and decreasing e 3 (see 
Fig. 8.1 in Ref. 


The ground state energy of Na + donors near the sur¬ 
face of Si was estimated in Ref. ISJ This case is com¬ 
plicated by an additional electric field in the inversion 
layer which tends to localize electrons near the interface 
Si/Si0 2 , to increase the donor ionization energy e\ and 
to decrease the decay length a of the wave function of a 
single donor. For experimental values of the electric field, 
the energy = 35 meV was obtained in Ref. [91 and [251 
Using Eq. Q and electron light mass mi ~ 0.2?n e we 
arrive at a = 23A. Again, due to the high experimental 
concentration n, this value can be considered to be a rea¬ 
sonable estimate from below for the experimental result 
a = 59A. 
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VI. CONCLUSION 

We have explored the surface conductivity of Si gated 
by an ionic liquid. We studied small induced concen¬ 
trations of holes n below 10 12 cm -2 and showed that at 
such concentrations the transport is due to the nearest 
neighbor hopping between acceptor like states formed by 
sparse excessive negative ions of the IL. This transport 
is similar to the nearest neighbor hopping between Na + 
donors in Na + implanted Si MOSFETs. The similarity 
is not only qualitative but also quantitative due to the 
similar average masses of holes and electrons in Si. This 
brings us to the conclusion that an ionic liquid acts as a 
set of sparse acceptors on the surface of Si. 

In the future it would be interesting to extend these 
results to smaller hole densities, where our theory should 
work better. One can study different ILs and see how 
the sizes of ions and dielectric constant of ILs affect hop¬ 
ping conductivity. One can also extend this work to new 
challenging materials using the same IL so that contrast¬ 


ing results with Si will help to understand better the 
parameters of these materials. For example, the exten¬ 
sion of measurements of surface conductivity in IL gated 
rubrene 3 i to lower hole concentrations may help to under¬ 
stand nature and parameters of charge carriers. 
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